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Abstract 

We present some straightforward applications of the QCD heavy quark expansion, 
stated in previous papers [1-3], to the inclusive widths of heavy flavour hadrons. We 
address the question of the Ds lifetime and argue that - barring Weak Annihilation 
(WA) - t{Ds) is expected to exceed t{D^) by several percent; on the other hand 
WA could provide a difference of up to 10 20% of any sign. We extract rric, mt 
and \Vcb\ from rgL(D+) and rsL(i?). The values of the quark masses are somewhat 
higher, but compatible with estimates from QCD sum rules; we obtain |V"cfe| — 0.043 
for t{B) = 1.4 psec and BRsl{B) = 10.5% . We discuss the associated uncertainties 
in the l/mq expansion as well as some consequences for other electroweak decays. 
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Our theoretical understanding of both exclusive and inclusive heavy flavour decays 
has improved considerably over the last few years. This progress has been driven 
largely by treatments that involve expanding weak transition amplitudes in terms of 
l/mg, niQ being the mass of the heavy flavour quark, i.e. mg = 'n^beauty; ''^charm- 
In this note we will analyze some phenomenological consequences that such a general 
treatment has for inclusive decays of beauty and charm hadrons. 

The paper will be organized as follows : in Sect.l we recapitulate briefly the salient 
features of our method which in Sect. 2 is then applied to the lifetime of Dg mesons; 
in Sect. 3 we extract the size of rric from the semileptonic D width and infer the value 
of rrih] the value of \Vcb\ is then determined from the semileptonic B width. In Sect. 4 
we present our conclusions. 

1 General Method 

On very general grounds one expects that the spectator ansatz correctly describes 
inclusive heavy flavour decays for asymptotically large values of rriQ; yet at finite mg 
there arise nonperturbative corrections to transition rates that among other things 
generate different lifetimes for the various types of hadrons Hq carrying the flavour 
Q. It was shown in refs. [1-5] how these preasymptotic effects can be incorporated 
in a self-consistent and systematic way; this will be exemplified now for the case of 
charm lifetimes. 

The width for the weak decay of a charmed hadron into an inclusive final state / 
is obtained from the transition operator that has been expanded into a series of local 
operators 

r = ^^\KM\'{coif)cc+^aaGc+J:,^{cr,q)iqr,c) + Oil/mt)} (1) 

where the dimensionless coefficients depend on the parton level characteristics 
of / and on the ratios of masses of the final state quarks to the mass of the c quark; 
KM denotes the appropriate product of the weak mixing angles. 

The operators that appear on the right hand side of eq.(l) represent a rather 
universal cast, e.g. they control semileptonic as well as nonleptonic transitions and 
also distributions like semileptonic spectra: 

• the scalar density cc describing the 'quasifree' decay of the c quark; 

• the chromomagnetic dipole operator ciaGc; 

• the four-fermion operator (cTq) (qTc) . It actually contains terms differing in both 
the flavour of the light quark q and the particular Lorentz structure of the 7-matrices 

The widths for Ds, D^, and Ac decays are then given by the diagonal matrix 
elements of the operator F from eq.(l) between the corresponding states Dg, D^, 
and Ac. These matrix elements depend on long distance bound state dynamics; thus 
they cannot be evaluated perturbatively. Yet this aspect of the strong interaction 

^This expansion is based on a sufficiently large energy release. It is not applicable for beauty decays into 
r leptons, whose partial width is to be calculated explicitly. It is, however, the relevant approach to other 
annihilation-induced semileptonic widths H. 
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cannot be ignored, for even the matrix element of the 'quasifree' operator cc contains 
nonperturbative effects, giving rise to some nontrivial width sphtting: 

{H,\cc\H,) = 2Mh^ + J-^{H,\acxGc\H,) - -L(^H,\c(tBfc\H,) + ^(l/m^) (2) 

with denoting the covariant derivative; a relativistic normahzation has been em- 
ployed here. In addition to the chromomagnetic operator another operator appears 
now, namely c(iD)^c which describes the kinetic energy of the charm quark in the 
gluon background field. 

The size of the matrix element of the chromomagnetic operator between the heavy 
flavour mesons can be extracted from the hyperfine spin splitting between the pseu- 
doscalar and vector states: 

^^{D\-czaGc\D) {Ml, - Ml) . (3) 

For baryons it vanishes: 

{A,\aaGc\Ac) ~ . (4) 

The following general conclusions can be drawn from this approach: 

• The leading nonperturbative corrections arise on the level. 

• There are two distinct sources of such corrections: (i) The matrix elements of 
the d = 5 operator ciaGc that appear in eq.(l); its weight depends on which particular 
channel is considered, (ii) The expansion of the matrix elements of the 'quasifree' 
operator cc, see eq.(2) that affects uniformly all decay modes of a particular hadron. 

• These corrections are flavour independent: they do not depend explicitly on the 
flavour of the spectator antiquark in the meson; they do, however, differentiate be- 
tween mesons and baryons. 

• Order I/ttLq corrections are produced by the four-fermion operators (crjg)(grjc) 
describing the nonspectator effects of 'Pauli Interference' (PI) and 'Weak Annihila- 
tion' (WA) or 'Weak Scattering' (WS) in baryons. These flavour-dependant operators 
generate different lifetimes even among the members of the same isomultipletQ. 

The analysis of I /ml [1-4] and 01 corrections shows that their size is 

in general quite large in charm decays. For example the chromomagnetic operator 
produces corrections ~ [3{Ml* — M|,)]/(2m^) ~ 0.4. Higher order corrections can 
thus be expected to be still significant and we have to be satisfied with a typically 
semiquantitative analysis. For example the theoretical predictions t{D^)/t{D^) ~ 2, 
BRgL(Z)"'") ~ 16% and BRgL(D°) ~ 8% have to be seen as in agreement with the 
experimental findings within the uncertainties of such a treatment. 

2 Ds vs. Lifetimes 

The lifetimes of and Dg mesons could a priori differ substantially from each 
other due to SU (3)pi breaking and in particular due to a different weight of WA in 

^Similar effects leading to the width splitting of the two eigenstates are described by the corresponding 
four fermion operator of the form (cu)^ . Its effect, however, is strongly suppressed and does not shift the 
central value of the width; we ignore it in our analysis. 
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r(D°) and T{Ds). A previous analysis [j^ suggested WA to be smallish already in 
charm meson decays; furthermore l/mg expansions naturally lead to the prediction 
of SU (3)pi breaking to be small in heavy flavour decays, namely of the order of several 
percent only Thus one predicts 

r(D,) r(D°) (5) 

to first approximation, in agreement with recent E687 data: 



r(DO) 



1.13 ±0.05 (6) 



A rather delicate analysis is required to go beyond the semi-quantitative prediction 
of eq.(5) to see whether, indeed, the Dg lifetime should be slightly longer than the 
lifetime, as suggested by present data, and by how much. 

There are two main sources for lifetime differences among charm mesons, namely 
explicitly flavour-dependent 1/ml terms and corrections of order '^s/ihadr/'^^c due to 
^[/(S)!?! breaking in the 'flavour- independent' 1/m^ contributions. 

Attempting to estimate relative corrections of a few percent in nonleptonic D 
decays is a bold enough enterprise and we will ignore processes contributing less 
than 1% to the total width. In that category are e.g. doubly Cabibbo suppressed 
transitions c — >■ dsu with a relative rate oc tan^ 6c — 3 - 10~^ as well as Penguin driven 
processes; the latter are suppressed by either tiny mixing angles with the b quark or 
by the small mass of the s quark. 

Numerically large effects can arise only from Cabibbo-allowed channels: WA is 
present in nonleptonic -D° decays, and in semileptonic as well as nonleptonic Dg 
decays; it also drives Dg tv. Cabibbo-suppressed modes c — > ssu, ddu can 
produce corrections to the total width of a few percent: WA affects -D°, and Dg 
decays; in addition to PI intervenes also in Dg modes. 

Since the impact of WA as compared to PI is reduced in meson decays @, 0, 0, 
it is natural to compare T{Dg) with t{D^) rather than with t{D~^). There are four 
distinct sources for a difference in T{Dg) vs. r{D^) exceeding the 1% level: 
(a) The decay Dg tu. (b) PI in those Cabibbo suppressed Dg decays that are 
driven by the quark level transition c — > ssd. (c) The effects of SU{3) breaking on 
the matrix elements of the chromomagnetic and kinetic energy operators, (d) WA in 
nonleptonic D^ and in nonleptonic as well as semileptonic Dg decays. 

While the corrections listed under (a) and (d) have been discussed extensively in 
the literature, the corrections referred to under (b) and especially (c) on the other 
hand are much less familiar. 

(a) The width for the decay Dg rvr is completely determined in terms of the 
axial decay constant for the Dg meson: 

nDg rv^) ^ ^"^""rPnMn. _ 

OTT 

For fi)^ ^ 210 MeV one gets numerically 

T{Dg ^ TUr) - .03 T{D^) . (8) 
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This effect necessarily reduces t{Ds) relative to t{D^). 

(b) PI in Dg appears in the c — > ssu channel. Its weight is expressed in terms of 
the matrix elements of the four-fermion operators 

{Bs\{cLlt,SL){sLl^CL)\Bs) , (5, | (cL7^y Sl) (sL7^y C^) (9) 

with known coefficients that are computed perturbatively, see refs.]^, the hybrid 
renormalization of these operators down from the scale has to be included. The 
most reliable way to estimate the effect of PI, we believe, is to relate it to the observed 
difference in the and widths. Both PI and WA contribute in reality to this 
difference; yet according to a detailed analysis 0, ^ |^ PI is the dominant effect in 
mesons (see also the discussion below). It is also worth noting that the size of the PI 
correction in as estimated theoretically reproduces the observed width splitting 
for reasonable values of fo - provided the hybrid renormalization of the operators 
involved is taken into account. 

It is easy to see that the structure of the operators responsible for PI in Dg is 
exactly the same as in if one replaces the d quark by the s quark and adds the 
extra factor tan^ 6c; it is then destructive as well. From the assumed dominance of PI 
in the D^-D^ lifetime difference one thus arrives at: 

5T;,^^{Ds) - S ■ tan2(^,)(r(D+) - T{D^)) ^ -S ■ 0.03r(D°) (10) 

where the factor S has been introduced to allow for SU (3) violation in the relevant 
matrix elements of the four-fermion operators 0. 

The factor S is expected to exceed unity somewhat; in the factorization approxima- 
tion it is given by the ratio (/d^/Zd)^. Various estimates yield the range S" = 1 -h 1.7; 
to be more specific we adopt S = 1.4. Then we conclude that 

5rint(/^s) ~ -o.o4r(D°) . (11) 



(c) As already stated in Sect.l the l/'m^ corrections are given by the appropriate 
expectation values of two dimension five operators. As far as the chromomagnetic 
operator is concerned one has the general expressions: 

^ (D^ldaGcl/}") ~ - ■ (M|,o* - Mlo) . (12a) 



2Mc ' ' 2 

3 



2Mr 



{Ds\ciaGc\Ds) - - ■ {Mfj, - M|, ) . (126) 



Since the measured values for D-D* and for the Dg-D* mass splittings are almost 
identical, the chromomagnetic operator cannot be expected to induce an appreciable 
difference between t{Ds) and t{D^). 

^SU{3)fi violation also affects the coefficients of the operators that depend on the mass of the quarks 
in the final state; those corrections, however, are of the order of {ms/nic)^ whereas the factor S is linear in 
ms//ijjg^(jj., and we neglect the former. 
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The observation that the hyperfine sphtting is largely independent of the flavour of 
the spectator can be understood in the following intuitive picture (see ref. 0): using 
a simple constituent description one flnds that the chromomagnetic fleld strength is 
proportional both to the chromomagnetic dipole moment of the spectator (or more 
generally, of the light degrees of freedom in the D meson) and to the wavefunction 
density at the origin, 1/ (r)^. The former is most naturally expected to decrease when 
the (current) mass of the spectator quark increases, whereas the latter is always as- 
sumed to increase when going from non-strange to strange particles. Such a behaviour 
can explicitly be demonstrated at least in the limit when the spectator becomes heavy 
enough as well 0. These two effects may thus offset each other. The conclusion about 
the equal strenght of the hyperflne splitting can and of course must be tested in B 
mesons where the mass formulae are more reliable due to the larger h quark mass. 

The second operator that generates corrections is the kinetic operator c(iD)^c 
which describes a nonrelativistic ( "Fermi" ) motion of the charm quark. As mentioned 
above one expects the spatial wavefunction to be more concentrated around the origin 
for Dg than for D mesons. This in turn implies that the mean value of is expected 
to be larger for Dg than for D mesons; in other words the charm quark undergoes 
more Fermi motion as a constituent of than of D mesons. Correspondingly the 
lifetime of the charm quark is prolonged by time dilation to a higher degree inside 
Dg than inside D mesons. Eq.(2) makes this connection quite explicit ^: the 
factor 1 — (p^)/2m^ appearing in the expression for the transition operator is actually 
nothing but the mean value of the Lorentz factor \/l — v'^ that suppresses the decay 
probablity of a particle in a moving frame. 

The qualitative trend of this effect is thus quite transparent. However its nu- 
merical size is not, at least not yet. One can expect future progress in QCD sum 
rules and/or simulations of QCD on the lattice to determine the appropriate matrix 
elements (ifc|c(^D)^c|ifc)- Yet for the problem at hand, namely the Dg-D^ lifetime 
difference one can in principle extract the relevant matrix element from the measured 
values of meson masses in the charm and beauty sector according to the following 



prescription |jTO|: 



^ {{Dg\c{iY)fc\Dg) - {D\c{iDfc\D)) 



2Md 
2mbm, 



-{[{MdJ - {Mn)] - [{Mb^ - {Mb)]] (13) 



rrib - rric 

where {Md^Ds,b,Bs) denote the 'spin averaged' meson masses: e.g. for D mesons 



(M) 



D 
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''The situation with the chromomagnetic dipole moment could in principle have been more complicated if 
the nonspectator (gluon) light degrees of freedom tended to form a ground state with a nonzero spin inside 
a meson; then for an arbitrarily weak interaction of the spectator's spin the chromomagnetic term would 
be still finite. This abstract possibility however contradicts the spectrum of states in heavy quarkonia - it 
would produce an additional hyperfine splitting pattern due to interaction of the heavy spins with the spin 
of the light degrees of freedom. 
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and likewise for the other mesons. Accordingly one finds 



A 10 MeV shift in any of the 'spin averaged' mass terms (M) in eq.(14) corresponds 
numerically to the kinetic energy operator generating approximately a 1% change in 
the ratio t{Ds)/t{D^). The meson masses have been measured with an accuracy of 2 
MeV or better, which is sufficient for our analysis - with the exception of the Bg-B* 



sector. Very recent LEP/CDF data |jTT| indicate that the SU{3) mass splitting in the 
beauty sector is practically the same as in the charm sector, namely Mb^ — M^o ~ 
(94.5±4.6MeV) vs. Md,-Md ^ (99.5d=0.7MeV). Nothing is known experimentally 
about the B* mass; on the other hand retaining only the leading non-trivial term in a 
heavy quark expansion, one would conclude that the hyperfine splitting in the Bg-B* 
system be the same as in the B^-B^ - simply because the analogous equality holds for 
charm mesons. Combining the preliminary measurement of the Bg mass with these 
theoretical expectations about the B* mass would lead to the rather surprising result 
that the mean momentum of the heavy quark is practically the same in strange and 
non-strange heavy-fiavour mesons; the Fermi motion of the charm quark could then 
cause a difference in t{Ds) vs. t(Z)°) of at most 1%. 

This conclusion can be confronted with a general expectation: it is natural to 
expect the value of (p^) to be at least of order of (400 Me V)^ even in non-strange 
heavy mesons]^; this by itself would suppress the width in charm by about 3% and 
increase the mass of the meson by 50 MeV. Assuming that SU{3) violation increases 
(p^) by a moderate factor 1.5 in the strange system one would then expect at least a 
2% decrease in the width of Dg. 

We think that this oversimplified yet transparent line of reasoning cannot be ruled 
out yet; it would seem quite premature to conclude that the Fermi motion plays a 
negligible role in the Dg-D^ lifetime difference. Namely the chromomagnetic field in 
Ds and D may not coinside so closely to predict the Bg mass with the necessary 
precision. For one has to allow for sizeable corrections to the mass formulae stated 
above where only the leading terms in l/rrib and l/rric have been retained since ttIc 
is not much larger than typical hadronic scales. There is no general reason to expect 
that the mass splitting formulae hold to better than 30% accuracy. Taking into 
account that the size of the hyperfine splitting in charm mesons constitutes about 
140 MeV, such an uncertainty translates into a 3-4% shift in the lifetime ratio, see 
eq.(14). Therefore we conclude: (a) The Fermi motion of the charm quark inside 
the meson may well prolong t{Ds) by typically a few percent relative to t{D). {j3) 
A better determination of the chromomagnetic field can be obtained from the B 
system. Therefore a more definite predictions can be made once the mass of the 
vector meson B* has been measured and the experimental uncertainty on M{Bs) has 
been decreased. 

One further comment is in order here. The discussion given above was concerned 
with nonleading corrections to the mass formulae for charm particles that could be 
sizeable because of the moderate value of rric. Could there be analogous corrections 

recent QCD sum rule estimate yielded about 0.5 ~ 0.6 GeV^ |jl2[. 
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to the expansion of the widths of charm particles? In particular there could be a 
sizeable deviation from the relation {Ds\ciaGc\Ds) — {D\ciaGc\D); this would show 
up as a violation of the prediction Mb* — Mb^ — Mb* — Mb that is suggested by 
a simple extrapolation from charm. This deviation would reflect the weight of non- 
leading corrections in the charm system and could a priori have a significant impact 
on the width; for the chromomagnetic interaction seems to be more important nu- 
merically than the kinetic energy However we do not expect this to amount to 
an important effect. For the chromomagnetic operator appears twice, namely in the 
expansion of the transition operator F, eq.(l), and in the 'quasi-free' operator cc; 
there is an almost complete cancellation between the chromomagnetic contributions 
from the nonleptonic and the two semileptonic channels if one uses for the coefficient 
C2 in eq.(l) the expression obtained in ref. Even allowing for a 30% uncertainty 
in the coefficient C2 the shift cannot exceed 1% for each 30 MeV in the hyperfine 
splitting of D or D^. (A more accurate prediction had to include the next to leading 
perturbative corrections in the coefficient C2; computing them represents a task that 
is straightforward in principle, although tedious in practice.) 

(d) The numerical impact of WA on charm meson lifetimes is the most obscure 
theoretical item in the analysis. The task becomes even harder for our present anal- 
ysis addressing a difference in the WA contribution to t{Ds) and to t{D^). The 
uncertainty centers mainly on the question of how much the WA amplitude suffers 
from helicity suppression. In the valence quark description the answer is easily given 
to lowest order: the WA rate is suppressed by the ratio m^/m^ where ruq denotes 
the largest quark mass in the final state. For a proper QCD treatment one has to 
use current rather than the larger constituent quark masses, at least for the 
corrections. That would mean that WA is negligible in decays where the appro- 
priate factor reads {ms/m^'^ ^ .01 and a fortiori in Dg decays where only non-strange 
quarks are present in the final state. (Semi-)IIard gluon radiation cannot circumvent 
this suppression ^ . For such gluon corrections - when properly accounted for - drive 
the hybrid renormalization of the corresponding four fermion operators which how- 
ever preserves their Lorentz structure and therefore does not eliminate the helicity 
suppression, at least in the leading log approximation. A helicity allowed amplitude 
can be induced only at the subleading as{ml)/n level and is thus expected to be 
numerically insignificant. 

On the other hand nonperturbative dynamics can quite naturally vitiate helicity 
suppression and thus provide the dominant source of WA. These nonperturbative 
effects enter through nonfactorizable contributions to the hadronic matrix elements. 
This has been analyzed in considerable technical detail in ref.[0. The expression 
for the weak annihilation operator in Dg as well as its renormalization can be easily 
obtained from the general expression for the case of or (see refs. 0, 0) by 
interchanging the colour factors Ci and C2: 



F 



Gp\Vcb 



{p- 



''^^^,u-p^,Pu){as^ng{cLl^,qL) {qHf^CL) + aoct{cL-f^,—qL) {qLlt^—CL)) ] 




ann — 



Gtt 



'smg 



(3c2 + 2ciC2)ae^/2 + \c\ - - l)(3c^ + 20,0^) 
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aoct = 2c?-(ae9/2_i)(2c^ + ^CiC2) , 

Ci ~ '^^ ' — ' Q = d for decays , 

ci ~ , C2 ~ , g = s for decays , 

^asiml)^ ' 3 ^ ^ ^ 

For each of the two semileptonic channels in Ds one additionally has 

«sing = ^a^'^' + I , floct = -^(«^'^' - 1) • (15a) 

Numerically in Dg the overall singlet coefficient is about 9 and the octet one is 
—2.3 whereas in they are —.32 and 3.2 respectively. 

The reason behind Oging being so much larger for Ds than for D^ mesons is that 
there is a colour- allowed WA contribution to Ds decays while WA is colour-suppressed 
in D^ decays. This colour- allowed contribution is obviously factorizable, however the 
factorizable part practically vanishes due to helicity suppression. Appreciable effects 
can then come only from nonfactorizable contributions or from OlasifriD/ir) cor- 
rections to the leading log approximation. Naive colour counting rules suggest that 
nonfactorizable parts in the matrix elements of colour singlet operators are sup- 
pressed as compared to those of colour octet operators. This line of reasoning is at 
best semi- quantitative, but if one adopts it one would conclude that the weight of 
WA is similar in inclusive Ds and D^ decays. As already stated nonleading pertur- 
bative corrections are capable of producing helicity unsuppressed contributions even 
to factorizable matrix elements; yet also they are colour-suppressed. 

It was shown in ref. that a detailed experimental study of the semileptonic 
width and the lepton spectrum, in particular in the endpoint region, in D° vs. D~^ 
and/or in vs. decays would allow us to extract size of the matrix elements 
that control the weight of WA in all inclusive B and D decays. Since such data are 
not (yet) available we can at present draw only a qualitative conclusion: WA is not 
expected to affect the total lifetimes of D^ and Ds mesons by more than 10 -r- 20%, 



see refs. P, Furthermore WA does not necessarily enhance t{Ds)/t{D^): due to 
its interference with the spectator reaction it could even reduce it! 

To summarize our findings on the Ds — D^ lifetime ratio: SU (3)pi breaking in 
the leading nonperturbative corrections of order can - due to 'time dilatation' 

- increase t{Ds) by 3 -i- 5%. On the l/m|? level there arise three additional effects. 
Destructive interference in Cabibbo suppressed Dg decays increases t{Ds) again by 
3 -7- 5% whereas the mode Dg tv decreases it by 3%. These three phenomena 
together lead to t{Ds)/t{D^) ~ 1.0 -r- 1.07. Any difference over and above that has 
to be attributed to WA. Taking these numbers at face value one can interprete the 



recent measurement ||14[ of the Ds lifetime in turn as more or less direct evidence 
for WA to contribute not more than 10-20 % of the lifetime ratio between charm 
mesons. As expected [|], 0, it does not constitute the major effect there. Finally the 
predictions just stated can be refined by future more accurate measurements, namely 
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• of the difference in the semileptonic spectra of charged and neutral mesons in the 
charm and in the beauty sector to extract the size of the matrix elements controlling 
the weight of WA; 

• of M(Ab), M{Bs) and M{B*) to better than 10 MeV to determine the expecta- 
tion values of the kinetic energy operator. 

3 Heavy Quark Masses and | Fc6 1 ~ phenomenological approach 

Conventional wisdom has it - based on considering the simplest perturbative diagrams 
- that semileptonic decays are easier to treat theoretically since they are less affected 
by the strong interactions. Our analysis of non-perturbative corrections in inclusive 
heavy flavour decays [1-5] offers additional evidence in support of this conviction. 
For there are smaller and fewer nonperturbative corrections in semileptonic than 
in nonleptonic decays of heavy flavour mesons. Even the first, 0{af) perturbative 
corrections may appear to be quite large in nonleptonic b decays owing to a relatively 
large mass of the charm quark in the final state Q. The semileptonic widths are then 
best suited to extract fundamental parameters like quark masses and KM angles. 
The semileptonic width of the D meson is given by the expansion 

^ = ^'^-1' ■ ^^''^^^ - ^^(-' - 25/4).«(x)). 

■(1 - Kd/ttiI + ^GD/ml) - z^{x)GD/ml + 0{al a J ml I /ml)} (16) 

where the phase space factors z account for the mass of the quark g = s, ci in the final 
state: 

2^o(a;) = 1 — 8x + — x"^ — 12x^ logx , zi{x) = {1 ~ x)^ , 

4'^(0) = 1 , 4'^(l) = 3/(27r2-25/2))~0.41 , x={mjm,)^ (17a) 

(the function Zq'^ can be found in ref. |]15|), while Kd and Gd denote the kinetic energy 
and the chromomagnetic matrix elements respectively: 

(176) 

It should be noted that the explicit form of the order perturbative correction in 
eq.(16) refers to the on shell (pole) definition for the charm quark mass. 

The expressions given above allow one to determine the mass of the c quark from 
the measured semileptonic width, provided the weak mixing angles \ Vcs\ and |Vcd| are 
known. We shall assume in our subsequent analysis that the weak mixing matrix is 
determined by the existence of only three generations; the quantities \ Vcs\ and \ Vcd\ are 
then known with an accuracy that is more than sufficient for our purposes]]. Adopting 

®We are grateful to V.Braun for presenting explicit arguments in favour of this option and to M.Shifman 
for a discussion of this problem. 

^There are some direct experimental measurements of these two angles, however they suffer from relatively 
large uncertainties for Vcs- 
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for the semileptonic widths of D mesons the value of FgL = BRgL(D+) ■ ~ 
1.06 ■ 10~^^ GeV and assuming (p^) ~ .3 GeV^ we then find 

mc = (1.57 ±0.03) GeV 

where we have included only the experimental error, coming mainly from FgL. We 
used here the values ^^(m^) = .33 and = 140 MeV. 

The smallness of the error in rric quoted above reflects the fact that the width 
depends on the fifth power of rric. The value found for rric is then also not very 
sensitive to details in the nonperturbative corrections to the spectator picture; these 
theoretical uncertainties will be discussed next. 

The perturbative 0{as) term reduces the width in eq.(16) by about 25%; the 
leading nonperturbative corrections ~ (9(l/m^) are of comparable size: the chromo- 
magnetic term and the kinetic term yield a reduction by ~ 25% and by about 6%, 
respectively. However all these effects do not generate a prominent change in the value 
of 771^. the perturbative corrections increase by 75 MeV; the chromomagnetic term 
and the kinetic energy term force up by 85 MeV and by 20 MeV, respectively, for 
the stated value of (p^). From these numbers we infer that the associated uncertain- 
ties in these corrections are rather insignificant. There is some uncertainty concerning 
the value of 05 (M|) and the scale at which to evaluate as in charm decays; yet those 
effects are quite unlikely to exceed 20% and can be addressed by including O^a^) 
contributions. Potentially larger errors can be expected from the nonperturbative 
effects. As discussed in the previous section the uncertainty in the chromomagnetic 
field could conceivably be of order 30% in the charm system. Corrections of similar 
size can be expected from higher power terms in the l/rric expansion. Finally the 
exact value of the kinetic term is not known. Yet its impact is generally somewhat 
suppressed as compared to the chromomagnetic interaction, and is dominated by the 



latter for reasonable sizes [|12| of the mean Fermi momentum of the heavy quark; the 
dependence on the (p^) is illustrated in Table 1. It is worth adding that for obvious 
reasons there is no significant dependence of the extracted value for rric on when 
the latter is varied within reasonable limits for a current quark mass. 

Combining all of this we then estimate the present theoretical uncertainty in ex- 
tracting rric to be about 30 MeV; to be conservative one may increase it up to say 
50 MeV: 

me = (1.57 ±0.03 ±0.05) GeV . (18) 

A more detailed understanding of the intrinsic limitations on the accuracy of such 
approach can be expected ^ in the future from explicit calculations of the higher 
order corrections, see e.g. refs. [|l^, |T7[] . 

Quite often another definition of the quark mass is used, namely the MS one 
corresponding to a Euclidean renormalization point = —toq. At the one- loop level 
they are related by 



,2 

3 tt' 

eq.(18) then yields 1.35 GeV for this definition of rric- 
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It might appear at this point that the fifth power dependence of the decay width on 
the mass of the quark allows one to make a rather accurate extraction of rric almost 
without any detailed information about the nature of nonperturbative corrections. 
Such a conclusion would however overstate the facts: for its validity rests on the 
absence of nonperturbative corrections of order l/mg in the total widths, as proven in 
refs. 0. This explained a posteriori why simple minded estimates made a long time 
ago that ignored nonperturbative corrections yielded a charm quark mass of around 
1.5 GeV. The real shortcoming of these estimates was the following more subtle 
point: in these models one cannot distinguish between the mass of the charm quark 
and of the charm hadron in an unambigous fashion; not surprisingly the estimates 
numerically fell somewhere in between. This alternative of course is a reformulation 
of the problem of nonperturbative l/rnq corrections to widths. The QCD approach 
ensures that if rric is understood as the (current) quark mass then these corrections 
are absent! This is a consequence of the conservation of the colour fiow in QCD as 



can be seen by simple quantum mechanical arguments (see ref. for details). 

Having extracted a value for the charm quark mass one can then determine the 
mass of the beauty quark by employing an expansion of the heavy fiavour hadron 



masses in terms of the heavy quark mass [10 



(20) 

unfortunately the accuracy of this expansion is obviously controlled by the charm 
quark mass. Note that the absence of the perturbative corrections on the left-hand 
side of eq.(20) implies that pole masses are assumed. For the same value of the 
kinetic term the central value appears to be 5.0 GeV and 4.5 GeV for the pole and the 
MS masses of the b quark, respectively. For such indirectly determined mass of the 
beauty quark the dependence on the size of (p^) is more significant as is illustrated 
by Table 1. 

The same consideration fixes also the value of the scale A that determines the 
asymptotic mass difference between the mass of the heavy flavour hadron and the 
mass of the constituent heavy quark. For the family of the lowest lying pseudoscalar 
and vector mesons one then has A ^ 300 MeV. This value is also correlated with the 
size of the kinetic term, and the latter seems to represent now the main uncertainty 



in the value of A. Note that we have shown in a recent paper [18| that there is no 



lower bound on A in contrast to a recent claim |jT9[; a priori A could have been even 
negative (see also |]17|). 

Above we have discussed only the most obvious uncertainties that one encounters 
in the numerical evaluation of the heavy quark mass. There is a number of additional 
purely theoretical corrections that in principle can affect them. Those corrections 



were discussed in some detail in a recent paper [0; a more comprehensive analysis 
will be given in a forthcoming publication |]T^. We therefore do not dwell on them 
here and only comment that they are not expected to affect significantly the values 
extracted for the masses. 

Turning to the extraction of the quark mixing parameter \Vcb\ we first note that 
phenomenological studies have shown that the total widths and lepton spectra in 
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semileptonic B decays depend mostly on the difference — rric rather than on the 
absolute values of the heavy quark masses. In the framework of the heavy quark 
expansion these findings are the reflection of the fact that the c quark is rather heavy 
even as seen from the scale of the b quark mass; they can be understood as meaning 
that the error one makes in using the measured masses of the charm and beauty 
hadrons rather than those of the charm and beauty quarks is reduced relative to its 
natural scale l/mq. This error can actually be reduced even further by using values 
for the quark masses as they were extracted above from the widths of D mesons. A 
value for \Vcb\ can then be obtained from the measured semileptonic width. 

From a theoretical point of view such an inclusive method has some clear advan- 
tages over an extraction of \Vcb\ from the exclusive decays B — »• D,D*lu suggested 
by the Heavy Quark Symmetry - even when taken at the "gold plated" point of zero 
recoil ||22| . For the accuracy of the symmetry is governed by the lightest quark mass 
in the process, namely on the other hand the expansion parameter for the 

inclusive widths is an inverse power of mf,, and this method can be applied even for 
light quarks in the final state. 

The necessary expression for the width is given by eq.(16) where the obvious 
substitutions h for c and c for s are made. Using the values tb — 1.4 ps and 
BRgL(i?) ^ .105 we arrive at the results shown in Table 1. The main depen- 
dence is again due to the kinetic term, however now it is rather weak; the uncertainty 
associated with the value of as is now smaller owing to the higher energy scale. It is 
tempting to conclude then that this method provides us at least three times better 
accuracy for the extraction of \Vcb\ than other methods discussed so far. 

The values for and that we have obtained in this Section are around the 
upper end of existing estimates. More conventional values emerge if one uses a smaller 
kinetic energy term of about (.4GeV)^. It is worth noting that the theoretical pre- 
dictions of the QCD sum rules for A prefer smaller values of about 400 MeV 
using this smaller value of A decreases the predicted value of the kinetic term 
On the other hand it is most natural to expect [|l^ the mean value of the Fermi mo- 
mentum not to exceed significantly A. Therefore we think that values for the h quark 
mass around 4.95 GeV together with a relatively small scale for the Fermi momenta 
of the order of 300 -r- 500 MeV are both selfconsistent and acceptable phenomenolog- 
ically. This hypothesis can and will be cross-checked in detailed study of spectra in 
semileptonic and radiative B decays. 
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4 Conclusion 

In this note we have applied a systematic expansion in 1/ mq that exists for the inclu- 
sive widths, to a few phenomenologically interesting issues concerning the properties 
of charm and beauty mesons. The general method we use has been suggested earlier 
||T], ^ 1^ and refined in subsequent papers ^, ^ |T^, [T^. The main object of our 



analysis was the pattern of the charm meson widths, including the Dg width. 

Generally the size of preasymptotic nonperturbative effects is almost of order unity 
in charm; this makes it difficult to arrive at conclusions that are both definite and 



'For the sake of definiteness we used rgj(6 ^ u) — .QlY'^'y{h 
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reliable. The detailed classification of the corrections enable us to conclude neverthe- 
less that the width splitting between Dg and must be reasonably small. There 
are four sources for this difference that can produce effects on the few percent level: 
1) A larger Lorentz reduction of the decay probability for the c quark inside the Dg 
meson due to a more rapid Fermi motion. 2) Destructive interference in the Cabibbo 
suppressed decays of Dg] each of these effects are estimated to decrease the width of 
Dg hy ^ -T- 5%. 3) The mode Dg — > TPr increases it by about 3%. 4) Nonleptonic 
decays of both mesons are affected by WA which is present also in semileptonic de- 
cays of Dg. According to the analysis of ref. ||^ all hadronic parameters governing 
the strength of WA in different decays can be obtained by a careful comparison of 
the lepton spectra in semileptonic decays of charged and neutral B mesons. However 
such data do not exist for now; therefore the size of WA remains an unknown and - 
unfortunately - the potentially largest effect numerically in t{Ds)/t{D^). 

Other lines of reasoning suggest that the impact of WA on inclusive charm decays 
can hardly be large: presumably it does not exceed 10 -r- 20% of the total width; 
furthermore it can be of either sign! We see therefore that the data on the lifetimes 
of charm mesons fit reasonably well the theoretically expected pattern. On the other 
hand the existing experimental determination of the Dg lifetime can be viewed as 
some phenomenological constraint on the non-factorizable matrix elements of the 
four-fermion operators in charm mesons: 

- + 0.034^1,1 - 0.39™! & 0,01 • (H22^)^ (21) 

Id Id 

where we have used the notations introduced in eq.(16) of ref. 0. Of course the 
extremely small number in the rhs of eq.(21) should not be taken too literally. 

The theoretical prediction for the Dg lifetime can be further clarified by accurate 
measurements of the masses of Bg and B* masses on the one hand and on the other 
a better understanding of the scale of Fermi motion in heavy hadrons which can be 
obtained from an accurate analysis of spectra in semileptonic and/or radiative decays 
of beauty particles. 

Similar preasymptotic corrections splitting the widths of beauty particles are ex- 
pected to be essentially smaller |^ . A careful analysis of the QCD corrections 0] 



has lead to the observation that the width difference between the two mass eigenstates 
in the Bg-Bg sector quite probably represents the largest numerical difference in the 
family of beauty mesons: 



T -■^^200MeV)2 ^^^^ 

This estimate is valid for fs^ - which acts as an expansion parameter - not too 
large. To a very good approximation one can identify the two mass eigenstates as 
CP eigenstates. Obviously the upper bound for the width difference is reached when 
all final states in the decay channel b —>■ ccs for the decays of Bg have the same CP 
parity and lAF^J ~ 2F(6 ccs). The rough estimate for this partial decay width 
is given by the parton expression and constitutes about 20%, therefore the estimate 
eq.(22) is sensible up to fs, ^ 300 MeV. 
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A straightforward, in principle, application of a l/mq expansion of inclusive decay 
widths [|1], 0] is the determination of the mass of charm quark from the semileptonic D 
width, then obtaining the b quark mass from the mass formulae and subsequently ex- 
tracting I V^fel from the B meson semileptonic width. In charm non-perturbative effects 
dominate the more familiar perturbative corrections: the former constitute about 35% 
whereas 0{as) corrections yield approximately 25%. None of these however produce 
a significant change in the resulting value of rric owing to the fifth power dependence 
of the decay width on the mass. In fact this conclusion is a consequence of the non- 
trivial fact that there are no corrections of order l/m^ to the heavy flavour widths, 
as shown in refs. [0, this is a peculiar feature of the gauge nature of the strong 
interaction producing the bound state in the initial state and driving hadronization 
dynamics in the final states, and it reflects the conservation of the colour current. 
A reasonable estimate for the theoretical uncertainty of the value of the c quark is 
about 50 MeV and one can count on it being decreased in the future. 

The masses of b and c quarks obtained in such a way naively seem to be some 
100 MeV higher than the conventional values inferred from the direct QCD analysis, 
and in particular from estimates based on charmonium sum rules. In fact a careful 
analysis^] undertaken recently suggested a rational explanation of this difference in 



the framework of QCD. We will report on it in the forthcoming paper |T7 . 

Using the observed semileptonic width of B mesons and including the leading 
corrections to the parton formulae one obtains the value of the KM mixing angle 
defining the strength of the 6 — > c transitions: 

1^61^ 0.043 . 

This method seems to be the most accurate and realiable way to obtain the value 
of \Vcb\- For its accuracy is governed by powers of I /nib as confronted to the decays 
into exclusive final state D and D* ||2^ where actually the c quark sets the mass 
scale for the corrections to the Heavy Spin-Flavour Symmetry. The corrections to the 
inclusive decays have been calculated explicitely through order 1 / and in principle 
the expansion can be extended, whereas effects for the exclusive decays 

may already constitute the limiting factor for improving the estimates. 

Our expressions for the semileptonic decay widths can be easily translated into 
the semi-phenomenological parameters Zc and that were introduced long ago to 
account for the final quark mass suppression in inclusive semileptonic widths as well 
as for all other possible corrections to the parton formulae. Their usual definition was 

rsL(&-g) = ^.-^i|^IK.P , g = u,c . (23) 

For the thus defined factor Zc we get values varying from .36 at small (p^) to .43 
for (p^) ~ 0.8 GeV^; these values seem to agree with the estimates obatined from 



experiment using fits based on some phenomenological models ||2^ of a heavy quark 
decay. The ratio Zu/ Zc — 2.08 appears to be independent of the scale of the Fermi 
motion, which also reproduces the expectations that one inferred from those models. 

^It was elaborated in joint discussions with M.Shifman and A.Vainshtein; we are grateful for their per- 
mission to mention it prior to publication of that result. 
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The phenomenological extraction of rric and mt, enables one to determine the 
hadronic parameter A which in the nonrelativistic description of a heavy hadron plays 
a role of the constituent mass of the spectator (s). As was stated in ref. A actually 
completely defines the leading, l/mg, non-perturbative shift in the average invariant 
mass of the final hadrons in semileptonic or radiative decays. In the quasi two particle 
decays like 6 ^ s + 7 the correction is given by Sm^p ~ Ami, + ^'(/i^g^jj.) if one 
neglects the mass of strange quark. We see that the average final state hadronic mass 
is increased by some 1.5 -i- 2 GeV^ as compared to a perturbative estimate. We note 
that once again the leading nonperturbative effects seem to be at least comparable 
in size to the perturbative corrections |^ even in beauty decays. Accordingly they 
have to be included in a proper quantitative treatment; together with the effect of 
the Fermi motion the A parameter describes the significant l/rrib corrections to the 
spectrum of photons, which in turn define the relative weight low lying exclusive final 
states can command in such decays. 



Note added: While this paper was written up we become aware ^Tj of the work of 
ref.||2^ which have a significant overlap with the part concerning the determination 
of quark masses and the mixing parameter. At that point we present a more detailed 
discussion of nonperturbative corrections as they are relevant for extracting rric, mb 
and |y(c6)|. The two treatments thus complement each other. As was mentioned in 
Section 3 we also believe that using the upper bound of Guralnik and Manohar - that 
constitutes an important element of the analysis of paper p8[ - is irrelevent 
this context. 
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(D|c(iD)2c|I3) 


rric, GeV 


rric (MS), GeV 


mb, GeV 


mfe(MS), GeV 


A, GeV 


\Vc,\ 


(0.1 GeV)^ 
(0.3 GeV)2 
(0.5 GeV)2 
(0.7GeV)2 
(0.8 GeV)2 


1.55 


1.33 


4.89 


4.43 


0.422 


0.0442 


1.56 


1.34 


4.91 


4.45 


0.390 


0.0437 


1.57 


1.35 


4.96 


4.50 


0.327 


0.0427 


1.59 


1.36 


5.03 


4.56 


0.234 


0.0413 


1.60 


1.37 


5.07 


4.60 


0.177 


0.0404 


(0.9 GeV)- 


1.(31 


1.39 


5.12 


1.61 


0.112 


0.0395 



Table 1: Dependence of the extracted parameters on the size of the kinetic energy 
operator in nonstrange mesons. We used here the strength of the QCD running 
couphng corresponding to Aqcj) = 180 MeV. 
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